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1. Introduction 

Nucleon electromagnetic form factors contain information about the size, shape and current 
distributions inside a nucleon. Conventionally the Dirac (F{) and Pauli (F 2 ) form factors are used 
to parametrize the nucleon electromagnetic matrix element through the following definition: 



{N(P'M m (x)\N(P)) 



y^F,{Q 2 ) + io^^-F 2 {Q 2 ) 



u(P), (1.1) 



where Q 2 = —q 2 = —{P' — P) 2 is the momentum transfer of the nucleon. Je, n (x) is the elec- 

'em,p 



tromagnetic current, with the explicit form j£ m „ = ^uy^u — kdy^d for the proton and Je mn 



— ^uj^u + ^dy^d for the neutron. 

Sachs electric (Ge) and magnetic (Gm) form factors are also frequently used by experimental- 
ists, and they are defined as the following linear combinations of the Dirac and Pauli form factors: 



Ce(Q : ) - F ]{ Q~)-^ I F 2 (Q 2 ) 1 (1.2) 
G M (Q 2 ) = F l (Q 2 )+F 2 (Q 2 ). (1.3) 



The isovector (Fj v 2 ) and isoscalar (F^ 2 ) form factors are defined, respectively, as 

K 2 (Q 2 ) = Ff^-F^Q 2 ) = Fl\ 2 {Q 2 )-Fi 2 (Q 2 )=F^ d (Q 2 ), (1.4) 
KiiQ 2 ) = Ff^+F^Q 2 ) = 1 (Fl 2 {Q 2 ) +Ft 2 (Q 2 )) = ^F"+ d (£2 2 ), (1.5) 

where F" 2 and Ff 2 are Dirac and Pauli form factors which parametrize the up and down quark 
vector currents, respectively. Calculations of the isoscalar form factors involve the evaluation of 
disconnected quark loops in the three-point correlation functions and are numerically expensive. 
In this talk we summarize the highlights of the recent results ^ for the isovector form factors, in 
which the disconnected quark loop contributions cancel in the exact isospin limit. In particular, 
we discuss the phenomenological dipole and tripole fits to the Q 2 dependence of the form factors, 
and study the chiral extrapolations using the formula from the small-scale expansion (SSE) and 
the covariant baryon chiral perturbation theory (BChPT). Similar studies with the mixed-action 
approach (domain wall valence on an Asqtad sea) can be found in Refs. [El 



2. Computational Details 

The calculations were performed on three 32 3 x 64 2+1 -flavor domain wall fermion gauge 
ensembles generated by the RBC and UKQCD collaborations [Q, |5|] with the input light quark 
masses set to ami = 0.004,0.006 and 0.008 and the strange quark mass fixed to am s = 0.03. The 
corresponding pion masses are roughly 297 MeV, 355 MeV and 403 MeV. The Iwasaki gauge 
action was used with /3 = 2.25, which gives a lattice spacing of a = 0.084 fm [§]. The extent of 
the fifth dimension, L s , was chosen to be 16, and the domain wall height was M5 = 1.8. The choice 
of these parameters gives rise to a residual mass of am res = 0.000660(8) in the chiral limit, which 
is about 1/6 of the lightest input quark mass. A coarse ensemble with j8 =2.13 and ami/am s = 
0.005/0.04 on the 24 3 x 64 lattice, with a pion mass of roughly 330 MeV, was also analyzed to 
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study the discretization errors. On this ensemble, the residual mass was determined to be am res = 
0.00315(1), and the lattice spacing was found to be a = 0.114 fm 

We computed the forward quark propagators with a Gaussian smeared source constructed 
from APE smeared gauge links. The parameters for the Gaussian smearing and the APE smearing 
were carefully tuned to minimize the overlap with the excited states and reduce the fluctuations 
from the source itself (see the Appendix in Ref. ^ for details.). On each gauge configuration 
we calculated forward quark propagators at time slices t = 10,26,42 and 58. The source-sink 
separation was chosen to be 12 for the fine ensembles, and 9 for the coarse ensemble, which both 
amount to about 1 fm in physical units. The coherent sink technique |Q] was employed to calculate 
the four sequential propagators associated with the four source locations simultaneously, leading 
to a factor of 4 reduction in the computation time. In this technique contaminations from other 
sinks are averaged to zero over the gauge configurations. We have verified that the results from 
the coherent-sink calculations were consistent with those from the independent-sink calculations, 
in which the sequential propagators were computed independently for each source. The form 
factors were obtained from the nucleon three-point functions using the overdetermined analysis 
as described in detail in Ref. [jl|]. And we obtained the vector current renormalization constant Zy 
by setting ZyF^O) = 1 for each ensemble. 

3. Phenomenological Fits to the Q 2 Dependence 

The Q 2 dependence of the isovector electromagnetic form factors is often found to be well 
described by a dipole form over a large Q 2 range. Since there is no theoretical foundation for such 
a Q 2 dependence, we performed fits to the Q 2 dependence of to the one-parameter dipole or 
tripole form of 

F?(Q 2 ) = - \ - 9N , n = 2,3, (3.1) 

{\ + q 2 /m 2 y 

while F 2 V was fit to the two-parameter dipole or tripole form of 

F 2 V (2 2 ) = t 2 \ , N , n = 2,3. (3.2) 

We saw no statistically significant difference between the dipole and tripole fits. We also investi- 
gated stability of the fits by varying the maximum Q 2 values included in the fits, as shown in Fig |[ 
While the fit quality decreases with larger Q 2 cutoffs, the fit parameters do not show statistically 
significant changes. Another feature of the dipole fits is that at small Q 2 values, the Dirac form 
factor appears to be lower than the dipole fit, while at large Q 2 , it tends to be higher, which is 
consistent with the observation of the experimental data [^, ||]. The Dirac and Pauli mean-squared 
radii, (r)') 2 and (r x 2 ) 2 , and the anomalous magnetic moment K v , were obtained from the dipole fits 
with a Q 2 cutoff of 0.5 GeV 2 through 

(rF) 2 = ^,i = l,2; k v = F1(0). (3.3) 

Similar behaviors of the Q 2 dependence were observed for the Sachs form factors G E {Q 2 ) and 
G^{Q 2 )- We show the results for G E (Q 2 ) from all the four ensembles in Fig. 2(a)| where the 



phenomenological parametrization [0] of the experimental results is also plotted for comparison. 
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Figure 1: Dipole fits to F{'(Q 2 ) and F{(Q 2 ) with three different cutoffs in Q 2 , 0.5, 0.7 and 1.1 GeV 2 , on the 
m K = 297 MeV ensemble. The top two graphs show the actual dipole fits to the lattice data. The bottom two 
graphs show the ratios of the data to the dipole fits with Q 2 < 0.5 GeV 2 . Fits with the other two fit ranges 
are normalized relative to the Q 2 < 0.5 GeV 2 fits. 



4. Chiral Extrapolations 

To compare the lattice results for the form factors at non-zero momentum transfer, we need to 
do chiral extrapolations at finite Q 2 values. Chiral perturbation theory requires that the momentum 
transfer values in the chiral expansions are small compared to the chiral scale of about 1 GeV. 
The available non-zero Q 2 values in our simulations range from 0.2 to 1.05 GeV 2 , which makes 
it unreliable to utilize the chiral formula for the Q 2 dependence of the form factors. Instead, we 
performed chiral fits to (r\ ) 2 , (rp 2 • k v (to get rid of the explicit k v dependence in the chiral formula) 
and K v . As discussed in Refs. [|l|, [Rj], the values of K v used in the chiral extrapolations should be 
rescaled with a factor of Mjf /M^ ys to get rid of the pion mass dependent magneton. We used two 
variants of the chiral effective field theories to perform the extrapolations. One is the heavy baryon 



chiral perturbation theory with the explicit A degrees of freedom Jl lQ, the so-called small-scale 
expansion (SSE), to <^(e 3 ). The other is the &{p A ) SU(2) covariant baryon chiral perturbation 
theory (BChPT) [12[ 13]. In order to disentangle the investigations of the applicability of the chiral 
effective field theories from the discretization effects, we included only the fine domain wall results 
in the chiral extrapolations discussed below. 

In the SSE chiral extrapolations, we performed simultaneous fits to {r\) 2 and (r^) 2 • K v (solid 
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Figure 2: (a) Comparison of lattice results for the isovector Sachs electric form factor G V E with the phe- 
nomenological parametrization of the experimental results (b) Chiral extrapolations of the nucleon mass 
using ff{p A ) SU(2) covariant baryon chiral perturbation theory. Only the a — 0.084 fm domain wall results 
(solid circles) are included in the fits. 



lines in Figs. 3(a) and 3(cJ), leaving the pion-nucleon-A coupling constant, ca, and a counter term, 
B r w (X), as free parameters, while other low-energy constants were fixed to their phenomenological 
values. Then the value of ca obtained from such fits was used as an input to determine three other 
unknown parameters (cy, K®,E[(X)) in K v (solid line in Fig. |3(e)[ ). As seen from the figures, the 
fits do not describe the lattice data very well. Adding a constant term to the formula for (r^) 2 • K v 
(dashed lines in Figs. 3(a) and 3(c) ) improves the fit quality, but the resulting extrapolated values 
at the physical point miss the empirical results [14, 15] by 10-20%. The difficulty is that the lattice 
data show a weaker curvature than the SSE expansion at the given order, which can be a result of 
i) the pion masses are still too heavy for the SSE formula to be accurate to the 10-20% level at the 
given order; or ii) the curvature of the data is obscured by un-controlled systematic errors, such as 
finite volume effects. 

A prerequisite to use the BChPT formula is to know the pion mass dependence of the nucleon 
mass, M^{m n ). We determined some of the low-energy constants needed in the &(p 4 ) BChPT 
formula for M/v(m^) from our three fine domain wall points, as shown in Fig. 2(b% In the BChPT 
chiral extrapolations, we performed simultaneous fits to {r\) 2 , (rJJ) 2 • K v and K v with four free pa- 
rameters. The resulting fits are shown in Figs. 3(bJ, 3(d)| and [3(f)| . Once again, the data show less 
curvature than the BChPT expansion at this given order. And the extrapolated physical values are 
also 10-20% lower than the empirical results. 



5. Conclusions 

We calculated the nucleon electromagnetic form factors with 2+1 flavors of domain wall 
fermions on three fine ensembles with a = 0.084 fm, and one coarse ensemble with a = 0.114 
fm. Focuses have been given to the study of the phenomenological dipole fits to the momentum 
transfer dependence of the isovector Dirac and Pauli form factors, and the investigations of chiral 
extrapolations to the isovector Dirac and Pauli mean-squared radii and the anomalous magnetic 
moment of the nucleon. We used two formulations of the baryon chiral effective field theories to 
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Figure 3: Chiral fits for the isovector Dirac and Pauli mean-squared radii and the anomalous magnetic 
moment of nucleon. The left panel shows the fits using the 0(e 3 ) SSE formula. The right panel shows the 
fits using the 0(p 4 ) BChPT formula. Only the a = 0.084 fm results are included in the fits. 



describe our data, the SSE formulation and the covariant baiyon chiral perturbation theory, and 
found that neither of these formulations can describe our data well. This may be caused by the 
relatively heavy pion masses in our simulations or the uncontrolled systematic errors, such as the 
finite volume effects at the lightest pion mass. To address these questions, we may need to simulate 
at several lighter pion masses at several different volumes to pin down the systematic errors. This 
remains a task for future investigations. 
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